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I explore the possibility that a quantum system S may be described completely by the combina¬ 
tion of its standard quantum state \ip) and a (hidden) quantum state \4>) (that lives in the same 
Hilbert space), such that the outcome of any standard projective measurement on the system S 
is determined once the two quantum states are specified. I construct an algorithm that retrieves 
the standard quantum-mechanical probabilities, which depend only on \tp), by assuming that the 
(hidden) quantum state \4>) is drawn at random from some fixed probability distribution Pr(.) and 
by averaging over Pr(.). Contextuality and Bell nonlocality turn out to emerge automatically from 
this algorithm as soon as the dimension of the Hilbert space of S is larger than 2. If \(f) is not com¬ 
pletely random, subtle testable deviations from standard quantum mechanics may arise in sequential 
measurements on single systems. 


I. INTRODUCTION 

Standard textbook quantum mechanics does not 
provide any explanation of what causes one particular 
outcome to occur in a given measurement on a quantum 
system S prepared in some state IV') G Hs- The quest 
for such an explanation is tantamount to a search for 
hidden variables [IH3- Whereas hidden variables are 
usually thought of as being classical, here I explore the 
possibility of using a hidden quantum variable, namely, 
a vector lying in the same Hilbert space Hs- (This 
Scune idea was first suggested by Bohm and Bub 
in 1966 ll^ , as Michael Hall pointed out to 
me. The algorithm 1 present is different from 
theirs, but the essence is the same.) I will 
describe an algorithm which assumes that the (hidden) 
vector is drawn at random from some fixed probability 
distribution Pr(.), such that the standard quantum state 
and the hidden quantum state together determine the 
outcome of any measurement one may perform, and such 
that the correct quantum-mechanical probabilities of 
the possible measurement outcomes are obtained when 
averaging over the probability distribution Pr(.). 

I apply the algorithm, first to measurements on finite¬ 
dimensional systems (a single qubit, a single qutrit, and a 
pair of qubits), and subsequently to infinite-dimensional 
systems (a radioactively decaying atom (or nucleus), and 
the two-slit experiment). 

II. ALGORITHM 

Here is an algorithm that takes as input the standard 
quantum state \tjj) of some system S, a hidden quantum 
state a probability distribution Pr(|^)) from which 
the hidden state has been drawn, and a standard projec¬ 
tive measurement M on S. The algorithm’s output con¬ 
sists of a single outcome for that measurement. Denote 
the Hilbert space dimension of system S by D, and as¬ 
sume that D is finite (the infinite-dimensional case is 
discussed below). In that case, there are exactly D out¬ 


comes for measurement M. 

1. Calculate the Born-rule probabilities for each of the 
outcomes of M from the state \4’) and sort them 
from largest to smallest. Write the sorted probabil¬ 
ities as = I in\ilf) p for n = 1... D. (That is, the 
eigenstates of M are now denoted by {|n)}„=i..._D-) 

2. Calculate corresponding Born-rule probabilities 
9n = I {n\4>o) P from the hidden state \(j)o)- 

3. Define tti such that the probability that qi is larger 
than TTi equals pi. In short-hand notation: 

Prob(gi > TTi) = pi. 

That is, define the subspace T of kets \(j)) such 
that I (lid)) P > TTi and adjust tti such that 

/j,Pr(|(()))d|(;)) =pi. 

4. Postulate that outcome 1 occurs iff qi > tti. 

5. Inductive rule: if outcome(s) 1... n did not oc¬ 
cur according to the algorithm, then postulate that 
outcome n -I- 1 occurs iff 

where 7r„+i is defined such that 

Prob'(g„+i > TTn+i) = - --, 

1 - Pi - P2 - - Pn 

where Prob’ denotes a conditional probability for 
qn+i, conditioned on the values of qi ...q^. The 
right-hand side is the conditional probability that 
outcome n-|-l occurs given that outcomes 1 through 
n did not occur. 

This algorithm, by construction, will produce outcome m 
with probability Pm = | P when one averages over 

Pr(.). Namely, by Rule 3, the probability with which 
outcome 1 is obtained equals pi. And by Rule 5, once we 
have shown that outcomes 1 through n occur with the 



correct probabilities pi.. .pn, then outcome n + 1 will 
occur with probability Pn+i- 

Two pure states are needed for having deterministic 
measurement outcomes. Whenever a system Si is entan¬ 
gled with some other system S 2 such that the joint state 
is pure but the state of Si is not, the algorithm is meant 
to be applied to the joint system. 

The two quantum states |'0) and |0o) do not play a 
symmetric role in the algorithm, because (i) the algo¬ 
rithm is constructed so as to reproduce the probabilities 
according to the standard quantum state \tjj) (but not 
those of |<('o)) and (ii) |(()o) is hidden and random. The 
standard quantum state is privileged in this sense. 

In the next Section I will illustrate the algorithm by 
examples that were constructed so as not to need the 
inductive Rule 5: it is always the most likely outcome 
(according to |' 0 )) that occurs. 


III. THE ALGORITHM AT WORK 

The algorithm does not specify the probability dis¬ 
tribution Pr(.) over the states \(j)) (it works for any 
distribution). For a Hilbert space of dimension 

D < 00 one natural choice is the uniform distribu¬ 
tion according to the unique unitarily invariant Haar 
measure. That is, the requirement on Pr(.) is that 
Pr(|^))d \(f)) = Pr([/ \(f)))d{U \(f))) for any unitary operator 
U acting on and any \(j}) € In this Section 

the uniform distribution is the only choice considered. 


A. One qubit 

The rule provided by the algorithm turns out to be 
simple for a 2D Hilbert space if we use the unitarily in¬ 
variant probability distribution Pr(|(/))). For any projec¬ 
tive measurement there are just two outcomes, and we 
denote the larger of the two corresponding probabilities 
by Pi = I (1|'0) P- The condition on 91 = | (I|^o) P for 
the corresponding outcome to occur is 

gi > TTi := 1 -pi. 

An alternative way of writing the condition for picking 
outcome 1 rather than 2, is -l-pi > q 2 +P 2 [lyl- This 
shows that the two quantum states do end up playing 
symmetric roles after all. As such, the two-dimensional 
case is similar to Spekkens’ toy model Q, in which an 
elementary physical system is modeled by assuming the 
underlying physical state is classical and contains two 
bits of information, but where one bit of information al¬ 
ways remains unknown. This restricted state of knowl¬ 
edge then turns out to mimic surprisingly many features 
of a quantum state. In the present model, on the other 
hand, there are two quantum bits of information, one of 
which is assumed unknown. This then reproduces all of 
the quantum mechanics of a qubit. 


B. Contextuality for a qutrit 

Let us move up one dimension and consider a qutrit. 
The most probable outcome of a measurement on a qutrit 
in a given state {-ip) occurs, according to the algorithm, 
when [ 2 ^ 

gi > TTi := I - y/pP. 

This condition is not symmetric between pi and qi . By 
way of example, consider a qutrit (viewed here as a spin-1 
particle, with the eigenstates of Sz chosen as basis states) 
in the state 

\iP) = V^\mz = -1)+V^\mz = 0)+y(L49|m^ = -f 1), 

and assume that the hidden state is 

{(p) = a/RS \mz = —l)-l--\/(L3 \mz = 0)-|--\/(L4 = -|-1). 

A measurement of Sz will yield the outcome ttIz = -I-I 
[2l|. On the other hand, if we measure a different observ¬ 
able, namely one with eigenstates |±) := (|toz = —I) ± 

I m, = 0))/V2,| TTiz = I), then the outcome will be |+) 
1^ . even though the outcome \mz = +1) is still a pos¬ 
sible outcome of this measurement, too. This feature 
is exactly the type of contextuality that appears in the 
Kochen-Specker theorem. That is, if this feature is ex¬ 
plicitly excluded as a possibility for a hidden variable the¬ 
ory, then the theorem shows that a hidden-variable model 
thus restricted cannot reproduce all quantum-mechanical 
predictions. 

C. Bell nonlocality for two qubits 

When considering joint measurements on two qubits in 
an entangled pure state they are to be treated as a single 
measurement on a 4D system. The condition on qi for 
outcome 1 to occur is now 

gi > TTi := 1 - ^/pP. 

Suppose we have two qubits, denoted by A and B, in the 
entangled state 

| 0 )a | 0 )B+yO | 1 )^ | 1 )b+v^ | 0 )^ 11)5 , 
and suppose that the hidden state is a product state 

= \^)a \^)B ’ 

with 

|<^)^ = 708 10)^ + 70211)^, 

and 

7)5 = 7010)5 + 70611 ) 3 . 

If we perform measurements on the pair of qubits in 
the standard basis (meaning that both A and B are 






measured in the basis {|0), |1)}) the outcome will be 
|0)^ |0 )b [23 ■ when we perform a different measure¬ 
ment on A, namely, one in the basis |±) = (|0) ± 11))/v^, 
while B is still measured in the standard basis, the joint 
measurement will have outcome |+)_ 4 | 1)5 so 

the measurement outcome on qubit B depends on what 
measurement is performed on A. This violates explicitly 
one of the assumptions that leads to Bell’s inequalities. 
The algorithm is thus explicitly nonlocal. (This very 
Scune observation was made about the Bohm-Bub 
model by Mattuck in 1981 [fl^ .) This nonlocality 
cannot be used for signaling as long as |(^o) is unknown. 
On the other hand, if one knew the hidden state |(/)o) 
one certainly could use that knowledge to achieve 
signaling. Hence, if Nature does not allow signaling, 
then a restriction must exist on how much knowledge of 
|(/)o) is possible. I leave it as an open question to find 
exactly what upper limit on knowledge of the hidden 
state is set by a no-signaling constraint. 


IV. THE ALGORITHM IN INFINITE 
DIMENSIONS 


The Haar measure does not exist in infinite¬ 
dimensional Hilbert spaces and we cannot pick a uni¬ 
form probability distribution over states in such a space. 
Instead, a fairly natural alternative can be given for a 
separable Hilbert space (with a countable basis): given 
a system S and its Hamiltonian, we could assign a prob¬ 
ability distribution over all possible states that depends 
explicitly on the energy, where, for example, we could as¬ 
sume uniform distributions (according to the correspond¬ 
ing Haar measures) on finite-dimensional degenerate sub¬ 
spaces of fixed energy. Alternatively, we could place an 
upper limit on the energy of the system and thereby de¬ 
fine a finite-dimensional Hilbert subspace. Then we could 
use the Haar measure on that subspace to define the uni¬ 
form distribution. The following examples are chosen so 
as to be largely independent of the choice of probability 
distribution for the hidden state. 

One sometimes uses inhnite-dimensional Hilbert spaces 
that are not separable. For example, the theoretical 
description of a spontaneously broken symmetry, as in 
the Higgs mechanism or the ferromagnetic phase of iron, 
does require the concept of unitarily inequivalent Hilbert 
spaces. The algorithm for such a case needs yet another 
technical adjustment [l^], but one still can define a prob¬ 
ability measure over all such inequivalent Hilbert spaces. 
And thus the algorithm presented here could attribute a 
spontaneously broken symmetry to an asymmetric hid¬ 
den quantum state. 


A. Phases of light and of Bose-Einstein 
condensates 

Consider an interference experiment in which two light 
beams with exactly N photons each impinge on a 50/50 
beamsplitter. Photodetectors (assumed perfect) placed 
at the two output ports of the beamsplitter count pho¬ 
tons. If in one output port N — n photons are counted 
then the other detector must count N + n photons (with 
—N < n < N). One may compare this experiment with 
a similar one with two light beams described by coherent 
states (“quasi-classical” states [ 13 ), both with amplitude 
|a| = y/N and a relative phase difference -d determined by 
sin i? = ^. In the latter case one expects unequal counts 
of photons in a ratio of approximately (n -|- N) l(n — N ), 
because the input beams have a well-defined phase differ¬ 
ence. On the other hand, in the experiment with num¬ 
ber states there is no such phase difference in the input 
state, but the measurement may very well produce un¬ 
equal counts nonetheless (see [l3 for a nice illustration 
how such a phase difference builds up during the mea¬ 
surement). 

With the algorithm this particular outcome would be 
explained by the hidden quantum state of two output 
light beams: for example, if that state (just before count¬ 
ing photons, but after the beamsplitter) were 

\(j)o) = 70)8 \N -n)\N + n) |0) |0) - |4) |2), 

(note there is no reason the hidden state should con¬ 
tain exactly 2N photons) then the explanation for the 
observed counting rates would be that the result was 
predetermined by the combination of the actual output 
state and the hidden state (especially its large first term). 
More generally, according to the algorithm, light beams 
may be said to contain a fixed number of photons and 
possess a well-defined phase difference (as revealed by an 
interference experiment), too. 

A similar story may be told about the phase of Bose- 
Einstein condensates [l3| . Even if two independently pre¬ 
pared condensates have no well-defined phase difference, 
according to the standard quantum description, an in¬ 
terference experiment will still display crisp interference 
fringes 0 as if a well-defined phase difference already 
existed. Again, according to the algorithm that measured 
phase difference is predetermined, thanks to the hidden 
state of the condensates. 

B. Radioactive decay Sz spontaneous emission 

A radioactively decaying nucleus or atom has always 
been the prime example of an indeterministic system, for 
which the quantum laws of nature seem to introduce ir¬ 
reducible and objective chances: The atom will decay 
at some point in time, quantum mechanics provides the 
probability distribution for those decay times, but noth¬ 
ing in the atom or anywhere else in the universe deter¬ 
mines what that point in time is going to be. 



The algorithm does say something about when the 
atom will decay. Since the atom itself will not remain 
in a pure state, one could assume monitoring not just 
the mass (or energy) of the atom, but the surroundings 
of the atom and detect the decay product. For example, 
in the case of spontaneous decay by emission of a pho¬ 
ton, the time of detection would be predetermined by the 
hidden state of the joint system consisting of the atom 
and the relevant modes of the electromagnetic field. The 
initial state of that joint system would be pure, with no 
photons present in the relevant modes and the atom in 
the excited state, and it would stay pme, as Mollow has 
shown in his analysis of this problem (l5j | . This pure state 
and the hidden state together determine when the atom 
decays. 


C. Two-slit experiment 

The mystery of the two-slit experiment lies in the com¬ 
bination of three propositions: (i) that single particles 
build up an interference pattern on a screen behind the 
two slits, (ii) that each particle must go through one slit 
only, but (hi) that with just a single slit open the inter¬ 
ference pattern disappears. 

Here is what the algorithm says about the two-slit ex¬ 
periment. Each particle will end up on a particular po¬ 
sition on the screen, as determined by its wave function 
'!/i(f^ and its hidden wave function 4>{f)\ repeating the ex¬ 
periment will produce the standard interference pattern, 
because although ' 0 (f) stays the same from one particle to 
the next, the hidden wave function varies randomly. The 
algorithm also gives a definite answer (either “slit 1 ” or 
“slit 2”) to the counterfactual question: “If we had mea¬ 
sured through which slit the particle went, what would 
we have found?” But, just as on the standard quantum- 
mechanical account, an actual measurement will destroy 
the interference pattern (because the wave function does 
collapse upon measurement). 

V. DISCUSSION 

Perhaps it is no surprise that one can concoct a hidden- 
variable model, necessarily nonlocal, contextual, and vi¬ 
olating the Pusey-Barrett-Rudolph criterion of indepen¬ 
dence, such that the hidden variables conspire to re¬ 
produce the probabilistic results of quantum mechanics. 
What is more surprising is that there exists a very sim¬ 
ple model achieving all this, in which a pair of quantum 
states—the standard state | 0 ) and a hidden state | 0 o)— 
fully describes the physical state of a system S in that 
the pair determines the outcomes of all measurements on 
S (including incompatible ones). 

The extension of standard quantum mechanics through 
the addition of a hidden quantum state evades no-go the¬ 
orems against particular types of hidden variables by vi¬ 
olating at least one of the underlying assumptions. As 


regards the three main no-go theorems, those by Kochen- 
Specker, by Bell, and by Pusey, Barrett and Rudolph, the 
algorithm has these three consequences/properties: (i) if 
the outcome of a measurement in a basis {|a), | 6 ), |c)} is 
I a), the measurement outcome in the basis {|a), |B), jC)} 
could be \B)] (ii) in a joint measurement on two systems 
A and B, the outcome of the measurement on A may 
depend on the measurement performed on system B; the 
algorithm is thus nonlocal; (hi) in a joint non-separable 
measurement on systems A and B, the hidden state may 
be entangled even if A and B were prepared indepen¬ 
dently. 

Bell inequality violations (as well as Hardy’s paradox 
[l^ and, by an obvious extension to measurements on 
three systems, the GHZ paradox [ 13 ) are all explained 
by the non-locality property (ii). 

A weak point of the algorithm as presented here is 
that it takes measurements as primitive, just as the stan¬ 
dard account of quantum mechanics does, too. On the 
other hand, compared to the many-worlds interpreta¬ 
tion, there is an economy of resources: there are just 
two non-branching wave functions as opposed to one gi¬ 
gantic branching wave function. Compared to Bohmian 
mechanics, which assign reality to both positions of cor¬ 
puscles and the wave function, the algorithm is more 
symmetric: it assigns reality to two wave functions. 

There are further avenues for exploration. I did not 
discuss how the hidden state evolves in time and to 
what state it collapses after a measurement (nor what 
precisely constitutes a measurement). One could assume 
it obeys the same Schrodinger equation as does |0) and 
that it collapses to a random state upon measurement 
(and that random state can be said to determine to 
what eigenstate |0) collapses). Alternatively, one could 
assume that the hidden quantum state is not completely 
random and not completely hidden. This could lead to 
subtle deviations from standard quantum-mechanical 
predictions, not concerning the probabilities of single 
measurement outcomes, but about sequences of mea¬ 
surements performed on a single system. For example, 
consider a sequence of measurements on a single spin-1/2 
system: first we measure Sz , then Sx , then Sz again, then 
Sx again., etcetera. Standard QM predicts that the mea¬ 
surement outcomes are truly random and uncorrelated. 
But if the hidden state |0o) is actually not completely 
randomized after each measurement, it might be that 
the Sx measurement outcomes are pairwise correlated 
(either positively or negatively), and similarly for the 
Sz measurement outcomes. (This observation, too, 
was made in a discussion of the Bohm-Bub model 
and an experiment was performed (in 1966!) to 
test this sort of prediction on photons (but 
not a single photon!) traveling through 3 
polarizers [ll 8 l] .) In any case, a no-signaling require¬ 
ment puts a restriction on how much information one 
may possess about \4>o)- 



Acknowledgements 


Bub model to me. 


I thank Tommaso Calarco for enjoyable and useful dis¬ 
cussions, and Michael Hall for pointing out the Bohm- 


[1] A. Einstein, B. Podolsky, and N. Rosen, Phys. Rev. 47, 
777 (1935). 

[2] J. S. Bell, Physics 1, 195 (1964). 

[3] J. S. Bell, Rev. Mod. Phys. 38, 447 (1966). 

[4] J. F. Clauser, M. A. Horne, A. Shimony, and R. A. Holt, 
Phys. Rev. Lett. 23, 880 (1969). 

[5] S. Kochen and E. Speaker, J. Math. Mech. 17, 59 (1967). 

[6] A. Peres, Journal of Physics A: Math. Gen. 24, L175 
(1991). 

[7] M. F. Pusey, J. Barrett, and T. Rudolph, Nature Physics 
8, 475 (2012). 

(81 D. Bohm and J. Bub, Reviews of Modern Physics 38, 
453 (1966). 

[9] R. W. Spekkens, Phys. Rev. A 75, 032110 (2007). 

[10] R. D. Mattuck, Physics Letters A 81, 331 (1981). 

[11] C. Cohen-Tannoudji, J. Dupont-Roc, and G. Grynberg, 
Photons and Atoms-Introduction to Quantum Electrody¬ 
namics, vol. 1 (1997). 

[12] K. Mplmer, Physical Review A 55, 3195 (1997). 

[13] J. Javanainen and S. M. Yoo, Phys. Rev. Lett. 76, 161 
(1996). 

[14] M. Andrews, C. Townsend, H.-J. Miesner, D. Durfee, 
D. Kurn, and W. Ketterle, Science 275, 637 (1997). 

[15] B. Mollow, Physical Review A 12, 1919 (1975). 

[16] L. Hardy, Phys. Rev. Lett. 68, 2981 (1992). 

[17] D. M. Greenberger, M. A. Horne, A. Shimony, and 
A. Zeilinger, Am. J. Phys 58, 1131 (1990). 


[18] C. Papaliolios, Physical Review Letters 18 , 622 (1967). 

[19] Here the set of states for which q\ -\-p\ = p 2 -\-q 2 is ignored 
as that set is of measure zero. Analogous sets of measure 
zero for the higher-dimensional cases are ignored, too. 

[20] For a D-dimensional Hilbert space tti = 1 — 

This follows from the fact that for a hxed ket jl) and 
for a ket \(f)) picked at random from the unitarily invari¬ 
ant probability distribution on a D-dimensional complex 
Hilbert space, Prob(| (1|<^) [^ > c) = (1 — c)^~^ for any 
constant c s.t. 0 < c < 1. 

[21] because pi = 0.49 and qi = 0.4 > 1 — = 0.3. 

[22] because pi « 0.51 and qi = 0.6 > 1 — « 0.29. 

[23] because pi = 0.5 and q\ = 0.32 > 1 — p\^^ « 0.21. 

[24] because pi — 0.45 and gi = 0.54 > 1 — p\^^ « 0.23. 

[25] For example, for the Heisenberg ferromagnet consisting 
of infinitely many spin-1/2 systems, we may parameterize 
the ground state manifold by two Euler angles, namely, 
by indicating the direction of each individual spin by 
those angles (and all spins in the ground state then point 
in that single direction). We could then postulate some 
distribution over the Euler angles. An excited state can 
be defined by flipping some finite number of spins of a 
given ground state. The hidden state can then be drawn 
from a distribution over ground states and the excited 
states thus defined. 




